Review of Big Bang Nucleosynthesis and 
Primordial Abundances 

David Tytler, John M. O'Meara, Nao Suzuki & Dan Lubin 
Center for Astrophysics and Space Sciences; 
University of California, San Diego; 
MS 0424; La Jolla; CA 92093-0424 

February 1, 2008 

Abstract 

Big Bang Nucleosynthesis (BBN) is the synthesis of the light nuclei, 
Deuterium, 3 He, 4 He and 7 Li during the first few minutes of the 
universe. This review concentrates on recent improvements in the 
measurement of the primordial (after BBN, and prior to modification) 
abundances of these nuclei. We mention improvement in the standard 
theory, and the non-standard extensions which are limited by the data. 

We have achieved an order of magnitude improvement in the pre- 
cision of the measurement of primordial D/H, using the HIRES spec- 
trograph on the W. M. Keck telescope to measure D in gas with very 
nearly primordial abundances towards quasars. From 1994 - 1996, it 
appeared that there could be a factor of ten range in primordial D/H, 
but today four examples of low D are secure. High D/H should be 
much easier to detect, and since there are no convincing examples, it 
must be extremely rare or non-existent. All data are consistent with 
a single low value for D/H, and the examples which are consistent 
with high D/H are readily interpreted as H contamination near the 
position of D. 

The new D/H measurements give the most accurate value for the 
baryon to photon ratio, r/, and hence the cosmological baryon density. 
A similar density is required to explain the amount of Lya absorption 
from neutral Hydrogen in the intergalactic medium (IGM) at redshift 



1 



z ~ 3, and to explain the fraction of baryons in local clusters of 
galaxies. 

The D/H measurements lead to predictions for the abundances of 
the other light nuclei, which generally agree with measurements. 

The remaining differences with some measurements can be ex- 
plained by a combination of measurement and analysis errors or changes 
in the abundances after BBN. The measurements do not require physics 
beyond the standard BBN model. Instead, the agreement between the 
abundances is used to limit the non-standard physics. 

New measurements are giving improved understanding of the diffi- 
culties in estimating the abundances of all the light nuclei, but unfor- 
tunately in most cases we are not yet seeing much improvement in the 
accuracy of the primordial abundances. Since we are now interested 
in the highest accuracy and reliability for all nuclei, the few objects 
with the most extensive observations give by far the most convincing 
results. 

Earlier measurements of 4 He may have obtained too low a value 
because the He emission line strengths were reduced by undetected 
stellar absorption lines. The systematic errors associated with the 4 He 
abundance have frequently been underestimated in the past, and this 
problem persists. When two groups use the same data and different 
ways to estimate the electron density and 4 He abundance, the results 
differ by more than the quoted systematic errors. While the methods 
used by Izotov & Thuan jjj seem to be an advance on those used 
before, the other method is reasonable, and hence the systematic error 
should encompass the range in results. 

The abundance of 7 Li is measured to high accuracy, but we do not 
know how much was produced prior to the formation of the stars, and 
how much was destroyed (depleted) in the stars. 6 Li helps limit the 
amount of depletion of 7 Li, but by an uncertain amount since it too 
has been depleted. 

BBN is successful because it uses known physics and measured 
cross-sections for the nuclear reactions. It gives accurate predictions 
for the abundances of five light nuclei as a function of the one free 
parameter r\. The other initial conditions seem natural: the universe 
began homogeneous and hotter than T > 10 11 K (30 Mev). The 
predicted abundances agree with most observations, and the required 
rj is consistent with other, less accurate, measurements of the baryon 
density. 
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New measurements of the baryon density, from the CMB, clus- 
ters of galaxies and the Lya forest, will give rj. Although the accu- 
racy might not exceed that obtained from D/H, this is an important 
advance because BBN then gives abundance predictions with no ad- 
justable parameters. 

New measurement in the coming years will give improved accuracy. 
Measurement of D/H in many more quasar spectra would improve 
the accuracy of D/H by a factor of a few, to a few percent, but even 
with improved methods of selecting the target quasars, this would 
need much more time on the largest telescopes. More reliable 4 He 
abundances might be obtained from spectra which have higher spectral 
and spatial resolution, to help correct for stellar absorption, higher 
signal to noise to show weaker emission lines, and more galaxies with 
low metal abundances, to minimize the extrapolation to primordial 
abundances. Measurements of 6 Li, Be and Boron in the same stars 
and observations of a variety of stars should give improved models for 
the depletion of 7 Li in halo stars, and hence tighter constraints on 
the primordial abundance. However, in general, it is hard to think of 
any new methods which could give any primordial abundances with 
an order of magnitude higher accuracy than those used today. This 
is a major unexploited opportunity, because it means that we can not 
yet test BBN to the accuracy of the predictions. 

1 Introduction 

There are now four main observations which validate the Big Bang the- 
ory: the expansion of the universe, the Planck spectrum of the Cosmic 
Microwave Background (CMB), the density fluctuations seen in the slight 
CMB anisotropy and in the local galaxy distribution, and BBN. Together, 
they show that the universe began hot and dense [|J. 

BBN occurs at the earliest times at which we have a detailed understand- 
ing of physical processes. It makes predictions which are relatively precise 
(10% - 0.1%), and which have been verified with a variety of data. It is crit- 
ically important that the standard theory (SBBN) predicts the abundances 
of several light nuclei (H, D, 3 He 4 He and 7 Li ) as a function of a single 
cosmological parameter, the baryon to photon ratio, r] = nb/n y [p| . The ratio 
of any two primordial abundances should give rj, and the measurement of the 
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other three tests the theory. 

The abundances of all the light elements have been measured in a number 
of terrestrial and astrophysical environments. Although it has often been 
hard to decide when these abundances are close to primordial, it has been 
clear for decades (e.g. [§J, ||) that there is general agreement with the BBN 
predictions for all the light nuclei. The main development in recent years 
has been the increased accuracy of measurement. In 1995 a factor of three 
range in the baryon density was considered Qb = 0.007 — 0.024. The low 
end of this range allowed no significant dark baryonic matter. Now the new 
D/H measurements towards quasars give = 0.019 ±0.0024 (95%) - a 13% 
error, and there have been improved measurements of the other nuclei. 



1.1 Other Reviews 



Many reviews of BBN have been published recently 
[Ol and 
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e.g. i 

some of which are lengthy: e.g. 



U, HI and 



It}. All modern cosmology texts contain a summary. Several recent books 



0, |1, 11 and 



contain the proceedings of meetings on this topic: 
The 1999 meeting of the International Astronomical Union (Symposium 198 
in Natal, Brazil) was on this topic, as are many reviews in upcoming special 
volumes of Physics Reports and New Astronomy, both in honor of the major 
contributions by David N. Schramm. 



2 Physics of BBN 



Excellent summaries are given in most books on cosmology e.g.: 
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23| , |24]] , and most of the reviews listed above, including [^J , and JTlJ . 



2.1 Historical Development 



The historical development of BBN is reviewed by flSB], [E7J , [H| and 
|| . Here we mention a few of the main events. 

The search for the origin of the elements lead to the modern Big Bang 
theory in the early 1950s. The expansion of the universe was widely accepted 
when Lemaitre J28| suggested that the universe began in an explosion of a 
dense unstable "primeval atom". By 1938 it was well established that the 
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abundances of the elements were similar in different astronomical locations, 
and hence potentially of cosmological significance. Gamow [29|, [30] asked 



whether nuclear reactions in the early universe might explain the abundances 
of the elements. This was the first examination of the physics of a dense 
expanding early universe, beyond the mathematical description of general 
relativity, and over the next few years this work developed into the modern 
big bang theory. Early models started with pure neutrons, and gave final 
abundances which depended on the unknown the density during BBN. Fermi 
& Turkevich showed that the lack of stable nuclei with mass 5 and 8 prevents 
significant production of nuclei more massive than 7 Li, leaving 4 He as the 
most abundant nucleus after H. Starting instead with all possible species, 
Hayashi |JT| first calculated the neutron to proton (n/p) ratio during BBN, 
and Alpher |52 realized that radiation would dominate the expansion. By 
1953 [33| the basic physics of BBN was in place. This work lead directly to 
the prediction of the CMB (e.g. Olive 1999b 0]), it explained the origin of 
D, and gave abundance predictions for 4 He similar to those obtained today 
with more accurate cross-sections. 

The predicted abundances have changed little in recent years, following 
earlier work by Peebles (1964) |9|, H °y le & Ta y ler ( 1964 ) E3> and Wagoner, 
Fowler & Hoyle (1967) f34|. The accuracy of the theory calculations have 
been improving, and they remain much more accurate than the measure- 
ments. For example, the fraction of the mass of all baryons which is 4 He, Y p , 



is predicted to within 5Y P < ±0.0002 [p5| . In a recent update, Buries et al. 



[H] uses Monte-Carlo realizations of reaction rates to find that the previous 
estimates of the uncertainties in the abundances for a given r\ were a factor 
of two too large. 



3 Key Physical Processes 
3.1 Baryogenesis 

The baryon to photon ratio rj is determined during baryogenesis 



It is not known when baryogenesis occurred. Sakharov [38] noted 
that three conditions are required: different interactions for matter and anti- 
matter (CP violation), interactions which change the baryon number, and 
departure from thermodynamic equilibrium. This last condition may be 
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satisfied in a first order phase transition, the GUT transition at 10 s, or 



perhaps the electroweak transition at 10 s. If baryogenesis occurred at 
the electroweak scale, then future measurements may lead to predictions for 
7], but if, alternatively, baryogenesis is at the GUT or inflation scale, it will 
be very hard to predict r\ (J. Ellis personal communication). 

The matter /anti- matter asymmetry of the universe (the r\ value) is at- 
tracting discussion in the popular science press because of the inauguration of 
major experiments to study CP violation in B mesons (|411, f42 ]; Economist, 
May 8 1999, 85-87). 



3.2 The main physical processes in BBN 

At early times, weak reactions keep the n/p ratio close to the equilibrium 
Boltzmann ratio. As the temperature, T, drops, n/p decreases. The n/p ratio 
is fixed ("frozen in") at a value of about 1/6 after the weak reaction rate is 
slower than the expansion rate. This is at about 1 second, when T ~ lMeV. 
The starting reaction n+p ~^_D + 7 makes D. At that time photodissociation 
of D is rapid because of the high entropy (low rj) and this prevents significant 
abundances of nuclei until, at 100 sec, the temperature has dropped to 0.1 
MeV, well below the binding energies of the light nuclei. About 20% of free 
neutrons decay prior to being incorporated into nuclei. The 4 He abundance 
is then given approximately by assuming that all remaining neutrons are 
incorporated into 4 He. 

The change in the abundances over time for one r] value is shown in Figure 
|l|, while the dependence of the final abundances on rj is shown in Figure 0, 
together with some recent measurements. 

In general, abundances are given by two cosmological parameters, the 
expansion rate and rj. Comparison with the strength of the weak reactions 
gives the n/p ratio, which determines Y p . Y p is relatively independent of rj 
because n/p depends on weak reactions between nucleons and leptons (not 
pairs of nucleons), and temperature. If rj is larger, nucleosynthesis starts 
earlier, more nucleons end up in 4 He, and Y p increases slightly. D and 3 He 
decrease simultaneously in compensation. Two channels contribute to the 
abundance of 7 Li in the rj range of interest, giving the same 7 Li for two 
values of rj. 
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4 Measurement of Primordial Abundances 



The goal is to measure the primordial abundance ratios of the light nuclei 
made in BBN. We normally measure the ratios of the abundances of two 
nuclei in the same gas, one of which is typically H, because it is the easiest 
to measure. 

The two main difficulties are the accuracy of the measurement and depar- 
tures from primordial abundances. The state of the art today (1 a) is about 
3% for Y p , 10% for D/H and 8% for 7 Li, for each object observed. These are 
random errors. The systematic errors are hard to estimate, usually unreli- 
able, and potentially much larger. 

By the earliest time at which we can observe objects, redshifts z ~ 6, we 
find heavy elements from stars in most gas. Although we expect that large 
volumes of the intergalactic medium (IGM) remain primordial today |43" |, we 
do not know how to obtain accurate abundances in this gas. Hence we must 
consider possible modifications of abundances. This is best done in gas with 
the lowest abundances of heavy elements, since this gas should have the least 
deviations caused by stars. 

The nuclei D, 3 He, 6 Li and 7 Li are all fragile and readily burned inside 
stars at relatively low temperatures of a few 10 6 K. They may appear depleted 
in the atmosphere of a star because the gas in the star has been above the 
critical temperature, and they will be depleted in the gas returned to the 
interstellar medium (ISM). Nuclei 3 He, 7 Li and especially 4 He are also made 
in stars. 



4.1 From Observed to Primordial Abundances 

Even when heavy element abundances are low, it is difficult to prove that 
abundances are primordial. Arguments include the following. 

Helium is observed in the ionized gas surrounding luminous young stars 
(H II regions), where O abundances are 0.02 to 0.2 times those in the sun. 
The 4 He mass fraction Y in different galaxies is plotted as a function of the 
abundance of O or N. The small change in Y with O or N is the clearest 
evidence that the Y is almost entirely primordial (e.g. Fig 2). Regression 
gives the predicted Y p for zero O or N [Q. The extrapolation is a small 
extension beyond the observed range, and the deduced primordial Y p is within 
the range of Y values for individual H II regions. The extrapolation should 
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be robust JIS , but some algorithms are sensitive to the few galaxies with the 
lowest metal abundances, which is dangerous because at least one of these 
values was underestimated by Olive, Skillman & Steigman ^ . 

For deuterium we use a similar argument. The observations are made 
in gas with two distinct metal abundances. The quasar absorbers have from 
0.01 to 0.001 of the solar C/H, while the ISM and pre-solar observations are 
near solar. Since D/H towards quasars is twice that in the ISM, 50% of the 
D is destroyed when abundances rise to near the solar level, and less than 
1% of D is expected to be destroyed in the quasar absorbers, much less than 
the random errors in individual measurements of D/H. Since there are no 
other known processes which destroy or make significant D (e.g.[|J, [47|), we 
should be observing primordial D/H in the quasar absorbers. 

Lithium is more problematic. Stars with a variety of low heavy element 
abundances (0.03 - 0.0003 of solar) show very similar abundances of 7 Li (PS"| 
Fig 3), which should be close to the primordial value. Some use the observed 
values in these "Spite plateau" stars as the BBN abundance, because of the 
small scatter and lack of variation with the abundances of other elements, 
but three factors should be considered. First, the detection of 6 Li in two of 
these stars suggests that both 6 Li and some 7 Li was been created prior to the 
formation of these stars. Second, the possible increase in the abundance of 
7 Li with the iron abundance also indicates that the 7 Li of the plateau stars 
is not primordial. If both the iron and the enhancement in the 7 Li have the 



same origin we could extrapolate back to zero metals f49], as for 4 He, but 
the enhanced 7 Li may come from cosmic ray interactions in the ISM, which 
makes extrapolation less reliable. Third, the amount of depletion is hard to 
estimate. Rotationally induced mixing has a small effect because there is 
little scatter on the Spite plateau, but other mechanisms may have depleted 
7 Li. In particular, gravitational settling should have occurred, and left less 
7 Li in the hotter plateau stars, but this is not seen, and we do not know why. 
More on this later. 

The primordial abundance of 3 He is the hardest to estimate, because 
stars are expected to both make and destroy this isotope, and there are no 
measurements in gas with abundances well below the solar value. 
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4.2 Key observational Requirements 



By way of introduction to the data, we list some of the key goals of ongoing 
measurements of the primordial abundances. 

• 4 He: High accuracy, robust measurement in a few places with the lowest 
metal abundances. 

• 3 He: Measurement in gas with much lower metal abundances, or an 
understanding of stellar production and destruction and the results of 
all stars integrated over the history of the Galaxy (Galactic chemical 
evolution). 

• D: The discovery of more quasar absorption systems with minimal H 
contamination. 

• 7 Li: Observations which determine the amount of depletion in halo 
stars, or which avoid this problem. Measurement of 6 Li, Be and B to 
help estimate production prior to halo star formation, and subsequent 
depletion. 

Since we are now obtaining "precision" measurements, it now seems best 
to make a few measurements with the highest possible accuracy and controls, 
in places with the least stellar processing, rather than multiple measurements 
of lower accuracy. We will now discuss observations of each of the nuclei, and 
especially D, in more detail. 



5 Deuterium in quasar spectra 

The abundance of deuterium (D or 2 H) is the most sensitive measure of the 
baryon density ||. No known processes make significant D, because it is so 
fragile (0, |p0| , and f52|). Gas ejected by stars should contain zero D, 
but substantial H, thus D/H decreases over time as more stars evolve and 
die. 

We can measure the primordial abundance in quasar spectra. The mea- 
surement is direct and accurate, and with one exception, simple. The com- 
plication is that the absorption by D is often contaminated or completely 
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obscured by the absorption from H, and even in the rare cases when contam- 
ination is small, superb spectra are required to distinguish D from H. 

Prior to the first detection of D in quasar spectra J55| , D/H was measured 
in the ISM and the solar system. The primordial abundance is larger, because 
D has been destroyed in stars. Though generally considered a factor of a few, 
some papers considered a factor of ten destruction |5^| . At that time, most 
measurements of 4 He gave low abundances, which predict a high primordial 
D/H, which would need to be depleted by a large factor to reach ISM values 
II- 

Reeves, Audouze, Fowler & Schramm [|J noted that the measurement 
of primordial D/H could provide an excellent estimate of the cosmological 
baryon density, and they used the ISM 3 He +D to concluded, with great 
caution, that primordial D/H was plausibly 7 + 3 xlO -5 . 

Adams suggested that it might be possible to measure primordial 
D/H towards low metallicity absorption line systems in the spectra of high 
redshift quasars. This gas is in the outer regions of galaxies or in the IGM, and 
it is not connected to the quasars. The importance of such measurements was 
well known in the field since late 1970s |57l| , but the task proved too difficult 



for 4-m class telescopes ([0, f59|, 0)- The high SNR QSO spectra obtained 
with the HIRES echelle spectrograph [|6l| on the W.M. Keck 10-m telescope 
provided the breakthrough. 

There are now three known absorption systems in which D/H is low: first, 
D/H = 3.24 ± 0.3 x 10~ 5 in the z a i, s = 3.572 Lyman limit absorption system 
(LLS) towards quasar 1937-1009 g|, @; second, D/H = 4.0±g;| x 10~ 5 



in the z a b s = 2.504 LLS towards quasar 1009+2956 |j3|, and third, D/H 



< 6.7 xlO towards quasar 0130-4021 [64]. This last case is the simplest 



found yet, and seems especially secure because the entire Lyman series is 
well fit by a single velocity component. The velocity of this component and 
its column density are well determined because many of its Lyman lines are 
unsaturated. Its Lya line is simple and symmetric, and can be fit using the 
H parameters determined by the other Lyman series lines, with no additional 
adjustments for the Lya absorption line. There is barely enough absorption 
at the expected position of D to allow low values of D/H, and there appears 
to be no possibility of high D/H. Indeed, the spectra of all three QSOs are 
inconsistent with high D/H. 

There remains uncertainty over a case at z a b s = 0.701 towards quasar 
1718+4807, because we lack spectra of the Lyman series lines which are 
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needed to determine the velocity distribution of the Hydrogen, and these 
spectra are of unusually low signal to noise, with about 200 times fewer 
photons per kms -1 than those from Keck. Webb et al. [B5|, [B6| assumed a 



single hydrogen component and found D/H = 25 ±5 x 10 5 , the best case for 
high D/H. Levshakov et al. |)7|] allow for non-Gaussian velocities and find 
D/H ~ 4.4 x 10" 5 , while Tytler et al. |§ find 8 x 10" 5 < D/H < 57 x 10' 



(95%) for a single Gaussian component, or D/H as low as zero if there are 
two hydrogen components, which is not unlikely. This quasar is then also 
consistent with low D/H. 

Recently Molaro et al. p5| claimed that D/H might be low in an absorber 
at z = 3.514 towards quasar APM 08279+5255, though they noted that 
higher D/H was also possible. Only one H I line, Lya, was used to estimate 
the hydrogen column density Nhi and we know that in such cases the column 
density can be highly uncertain. Their Figure 1 (panels a and b) shows that 
there is a tiny difference between D/H = 1.5 xl0~ 5 and 21 xlO -5 , and it is 
clear that much lower D is also acceptable because there can be H additional 
contamination in the D region of the spectrum. Levshakov et al. [^] show 
that Nhi= 15.7 (too low to show D) gives an excellent fit to these spectra, and 
they argue that this is a more realistic result because the metal abundances 
and temperatures are then normal, rather than being anomalously low with 
the high Nhi preferred by Molaro et al. 

The first to publish a D/H estimate using high signal to noise spectra 
from the Keck telescope with the HIRES spectrograph were Songaila et al. 
[[HH , who reported an upper limit of D/H < 25 x 10~ 5 in the z ahs = 3.32 
Lyman limit system (LLS) towards quasar 0014+813. Using different spectra, 
Carswell et al. |6(| reported < 60 x 10~ 5 in the same object, and they found 
no reason to think that the deuterium abundance might be as high as their 
limit. Improved spectra ]7^j support the early conclusions: D/H < 35 x 10 -5 
for this quasar. High D/H is allowed, but is highly unlikely because the 
absorption near D is at the wrong velocity, by 17 ± 2 km s" 1 , it is too wide, 
and it does not have the expected distribution of absorption in velocity, which 
is given by the H absorption. Instead this absorption is readily explained 
entirely by H (D/H ~ 0) at a different redshift. 

Very few LLS have a velocity structure simple enough to show deuterium. 
Absorption by H usually absorbs most of the quasar flux near where the D 
line is expected, and hence we obtain no information of the column density 
of D. In these extremely common cases, very high D/H is allowed, but only 
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because we have essentially no information. 

All quasar spectra are consistent with low primordial D/H ratio, D/H ~ 
3.4 x 10" 5 . Two quasars (1937-1009 & 1009+2956) are inconsistent with D/H 
> 5 x 10" 5 , and the third (0130-4021) is inconsistent with D/H > 6.7 x 10" 5 . 
Hence D/H is low in these three places. Several quasars allow high D/H, but 
in all cases this can be explained by contamination by H, which we discuss 
more below, because this is the key topic of controversy. 



5.1 ISM D/H 



Observations of D in the ISM are reviewed by Lemoine et al. |73|. The first 
measurement in the ISM, D/H = 1.4 ± 0.2 xl0~ 5 , using Lyman absorption 
lines observed with the Copernicus satellite f74|, have been confirmed with 
superior HST spectra. A major program by Linsky et al 
a secure value for local ISM (< 20 pc) D/H = 1.6 ± 0.1 xlO 



75| , |76fl has given 

5 



Some measurements have indicated variation, and especially low D/H, 

71. Vidal- 



in the local and more distant ISM towards a few stars 
Madjar & Gry [^] concluded that the different lines of sight gave different 
D/H, but those early data may have been inadequate to quantify complex 
velocity structure |f7|]. Variation is expected, but at a low level, from different 
amounts of stellar processing and infall of IGM gas, which leaves differing 
D/H if the gas is not mixed in a large volume. 

Lemoine et al. [78 suggested variation of D/H towards G191-B2B, while 
Vidal-Madjar et al. 
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described the variation as real, however new STIS 
spectra do not confirm this, and give the usual D/H value. The STIS spectra 
jSOU show a simpler velocity structure, and a lower flux at the D velocity, 
perhaps because of difficulties with the background subtraction in the GHRS 
spectra. 

Hebrand et al. |8T| report the possibility of low D/H < 1.6 x 10~ 5 towards 
Sirius A, B. 

The only other instance of unusually low D/H from recent data is D/H 
5 (90%) towards the star 5 Ori |E2 



74+ - 19 

u - '^-0.13 



xlO" 



We would much like 
to see improved data on this star, because a new instrument was used, the 
signal to noise is very low, and the velocity distribution of the D had to be 
taken from the N I line, rather than from the H I. 

Possible variations in D/H in the local ISM have no obvious connections 
to the D/H towards quasars, where the absorbing clouds are 100 times larger, 
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in the outer halos young of galaxies rather in the dense disk, and the influence 
of stars should be slight because heavy element abundances are 100 to 1000 
times smaller. 

Chengalur, Braun & Burton |83j report D/H = 3.9 ± 1.0 xlO -5 from 
the marginal detection of radio emission from the hyper-fine transition of D 
at 327 MHz (92 cm). This observation was of the ISM in the direction of 
the Galactic anti-center, where the molecular column density is low, so that 
most D should be atomic. The D/H is higher than in the local ISM, and 
similar to the primordial value, as expected, because there has been little 
stellar processing in this direction. 

Deuterium has been detected in molecules in the ISM. Some of these 
results are considered less secure because of fractionation and in low density 
regions, HD is more readily destroyed by ultraviolet radiation, because its 
abundance is too low to provide self shielding, making HD/H 2 smaller than 
D/H. 

However, Wright et al. |84| deduce D/H = 1.0 ±0.3 x 10~ 5 from the first 



detection of the 112 /mi pure rotation line of HD outside the solar system, 
towards the dense warm molecular clouds in the Orion bar, where most D 
is expected to be in HD, so that D/H ~ HD/H2. This D/H is low, but not 
significantly lower than in the local ISM, especially because the H2 column 
density was hard to measure. 

Lubowichet al. [g|, H report D/H = 0.2±1 xl0~ 5 from DCN in the Sgr 



A molecular cloud near the Galactic center, later revised to 0.3 x 10~ 5 (private 
communication 1999). This detection has two important implications. First, 
there must be a source of D, because all of the gas here should have been 
inside at least one star, leaving no detectable D. Nucleosynthesis is ruled out 
because this would enhance the Li and B abundances by orders of magnitude, 
contrary to observations. Infall of less processed gas seems likely. Second, 
the low D/H in the Galactic center implies that there is no major source 
of D, otherwise D/H could be very high. However, this is not completely 
secure, since we could imagine a fortuitous cancellation between creation 
and destruction of D. 

We eagerly anticipate a dramatic improvement in the data on the ISM 
in the coming years. The FUSE satellite, launched in 1999, will measure the 
D and H Lyman lines towards thousands of stars and a few quasars, while 
SOFIA (2002) and FIRST (2007) will measure HD in dense molecular clouds. 
The new GMAT radio telescope should allow secure detection of D 82 cm 
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emission from the outer Galaxy, while the Square Kilometer Array Interfer- 
ometer would be able to image this D emission in the outer regions of nearby 
galaxies; regions with low metal abundances. These data should give the re- 
lationship between metal abundance and D/H, and especially determine the 
fluctuations of D/H at a given metal abundance which will better determine 
Galactic chemical evolution, and, we expect, allow an accurate prediction of 
primordial D/H independent of the QSO observations. 



5.2 Solar System D/H 

The D/H in the ISM from which the solar system formed 4.6 Gyr ago can 
be deduced from the D in the solar system today, since there should be no 
change in D/H, except in the sun. 

Measurement in the atmosphere of Jupiter will give the pre-solar D/H pro- 
vided (1) most of Jupiter's mass was accreted directly from the gas phase, and 
not from icy planetessimals, which, like comets today, have excess D/H by 
fractionation, and (2) the unknown mechanisms which deplete He in Jupiter's 
atmosphere do not depend on mass. Mahaffy et al. [87] find D/H = 2.6±0.7 



xl0~ 5 from the Galileo probe mass spectrometer. Feuchtgruber et al. 
used infrared spectra of the pure rotational lines of HD at 37.7 /im to measure 
D/H = S.SiiJxIO -5 in Uranus and 6.5ii;5xlCT 5 in Neptune, which are both 
sensibly higher because these planets are known to be primarily composed 
of ices which have excess D/H. 

The pre-solar D/H can also be deduced indirectly from the present solar 
wind, assuming that the pre-solar D was converted into 3 He. The present 
3 He/ 4 He ratio is measured and corrected for (1) changes in 3 He/H and 
4 He/H because of burning in the sun, (2) the changes in isotope ratios in the 
chromosphere and corona, and (3) the 3 He present in the pre-solar gas. Geiss 
& Gloeckler |89| reported D/H = 2.1±0.5 x!0~ 5 , later revised to 1.94±0.36 



xl(T 5 [g. 

The present ISM D/H = 1.6 ± 0.1 xlCT 5 is lower, as expected, and con- 
sistent with Galactic chemical evolution models, which we now mention. 



5.3 Galactic Chemical Evolution of D 

Numerical models are constructed to follow the evolution of the abundances 
of the elements in the ISM of our Galaxy. 
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The main parameters of the model include the yields of different stars, 
the distribution of stellar masses, the star formation rate, and the infall and 
outflow of gas. These parameters are adjusted to fit many different data. 
These Galactic chemical evolution models are especially useful to compare 
abundances at different epochs, for example, D/H today, in the ISM when 
the solar system formed, and primordially. 



In an analysis of a variety of different models, Tosi et al. |)Tfl concluded 
that the destruction of D in our Galaxy was at most a factor of a few, consis- 
tent with low but not high primordial D. They find that all models, which are 
consistent with all Galactic data, destroy D in the ISM today by less than a 
factor of three. Such chemical evolution will destroy an insignificant amount 
of D when metal abundances low as seen in the quasar absorbers. 



Others have designed models which do destroy more D , |K| , |J4 
for example, by cycling most gas through low mass stars and removing the 
metals made by the accompanying high mass stars from the Galaxy. These 
models were designed to reduce high primordial D/H, expected from the low 
Y p values prevalent at that time, to the low ISM values. Tosi et al. |)T 



describe the generic difficulties with these models. To destroy 90% of the D, 
90% of the gas must have been processed in and ejected from stars. These 
stars would then release more metals than are seen. If the gas is removed 
(e.g. expelled from the galaxy) to hide the metals, then the ratio of the 
mass in gas to that in remnants is would be lower than observed. Infall of 
primordial gas does not help, because this brings in excess D. These models 
also fail to deplete the D in quasar absorbers, because the stars which deplete 
the D, by ejecting gas without D, also eject carbon. The low abundance of 
carbon in the absorbers limits the destruction of D to <1% 



5.4 Questions About D/H 

Here we review some common questions about D/H in quasar spectra. 



5.4.1 Why is saturation of absorption lines important? 

Wampler |)5|] suggested that the low D/H values might be inaccurate because 
in some cases the H absorption lines have zero flux in their cores; they are 
saturated. Songaila, Wampler & Cowie |9(J suggested that this well known 
problem might lead to errors in the H column density, but later work, using 
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better data and more detailed analyses |9?[ has shown that these concerns 
were not significant, and that the initial result [98] was reliable. 

Neutral deuterium (D I) is detected in Lyman series absorption lines, 
which are adjacent to the H I lines. The separation of 82 km s _1 is easily 
resolved in high resolution spectra, but it is not enough to move D out of the 
absorption by the H. The Lyman series lines lie between 1216A and 9 12 A, 
and can be observed from the ground at redshifts > 2.5. 

Ideally, many (in the best cases > 20) Lyman lines are observed, to help 
determine the column density (Nhi, measured in H I atoms per cm -2 along 
the line of sight) and velocity width (b values, b = \/2cr, measured in km s _1 ) 
of the H. But in some cases only Lya has been observed (Q1718+4807, APM 
08279+5255), and these give highly uncertain D/H, or no useful information. 

The column densities of H and D are estimated from the precise shapes 
of their absorption lines in the spectra. For H, the main difficulties are the 
accuracy of the column density and the measurement of the distribution in 
velocity of this H. For D the main problem is contamination by H, which we 
discuss below. 

It is well known that column densities are harder to measure when ab- 
sorption lines become saturated. The amount of absorption increases linearly 
with the column density as long as only a small fraction of the photons at 
the line central wavelength are absorbed. Lines saturate when most photons 
are absorbed. The amount of absorption then increases with the log of the 
column density. 

Wampler |95] has suggested that D/H values could be 3 - 4 times higher 
in Q1937-1009 than measured by Tytler, Fan & Buries [53|. He argued that 
saturation of the H Lyman series lines could allow lower Nhi- This would 
lead to residual flux in the Lyman continuum, which would contradict the 
data, but Wampler suggested that the background subtraction might have 
been faulty, which was not a known problem with HIRES. 

Tytler Sz Buries explained why Wampler's general concerns were not 
applicable to the existing data on Q1937-1009. Thirteen Lyman series lines 
were observed and used to obtain the Nhi- The cross section for absorption 
(oscillator strength) decreases by 2000 from the Lya to the Ly-19 line. This 
means that the lines vary significantly in shape, and this is readily seen 
in spectra with high resolution and high signal to noise. The background 
subtraction looked excellent because the line cores were near zero flux, as 
expected. 
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Songaila, Wampler & Cowie |96|] measured the residual flux in the Lyman 
continuum of the D/H absorber in Q1937-1009. They found a lower Nhi and 
hence a higher D/H. Buries & Tytler |J7[] presented a more detailed analysis 
of better data, and found a lower Nhi, consistent with that obtained from 
the fitting of Lyman series lines. They explained that Songaila, Wampler 
and Cowie |96| had underestimated Nhi because they used poor estimates 
of the continuum level and the flux in the Lyman continuum. 

In summary, saturation does make the estimation of Nhi harder. Column 
densities of H might be unreliable in data with low spectral resolution, or low 
signal to noise, and when only a few Lyman lines are observed. The above 
studies show that it is not a problem with the data available on Q1937-1009, 
Q1009+2956, Q0014+8118 and Q0130-4021. For the first two quasars, we 
obtain the same answer by two independent methods, and for the last three 
the higher order Lyman lines are not saturated. 

Saturation is avoided in absorbers with lower Nhi, but then the D lines 
are weaker, and contamination by H lines becomes the dominant problem. 



5.4.2 Hidden Velocity Structure 

To obtain D/H we need to estimate the column densities of D and H. Column 
densities depend on velocity distributions, and when lines are saturated, it 
is hard to deduce these velocity distributions. Similar line profiles are made 
when the velocity dispersion is increased to compensate for a decrease in the 
column density. We mentioned above that this degeneracy is broken when we 
observe lines along the Lyman series. For Q1937-1009, which has the most 
saturated H lines of the quasars under discussion, Buries &; Tytler |^] showed 
that the D/H did not change for arbitrary velocity structures, constrained 
only by the spectra. The same conclusion was obtained for Q1009+2956 
[ |63|j . The favorable results for these two quasars do not mean that we will 
always be able to break the degeneracy. That must be determined for each 
absorption system. 

There are two reasons why hidden velocity structure is not expected to be 
a major problem. First, we are concerned about hidden components which 
have high columns and low enough velocity dispersions that they hide inside 
the wider lines from lower column gas. Such gas would be seen in other 
lines which are not saturated: the D lines and the metal lines from ions 
with similar (low) ionization. Second, we search for D in absorbers with 
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the simplest velocity distributions. They tend to have both narrow overall 
velocity widths and low temperatures, which makes it much harder to hide 
unseen components. Typically, the main component accounts for all of the 
absorption in the higher order Lyman lines, and these lines are too narrow 
for significant hidden absorption. 



5.5 Correlated Velocity Structure: Mesoturbulence 



In a series of papers, Levshakov, Kegel & Takahara p9[ , ||100|| , ||101|| have 
demonstrated a viable alternative model for the velocity distribution. 

In most papers, absorption lines are modelled by Voigt profiles. The line 
width is the sum of the thermal broadening, turbulent broadening, and the 
instrumental resolution, each of which is assumed to be Gaussian. When an 
absorption line is more complex than a single Voigt, gas centered at other 
velocities is added to the model. As the signal to noise increases, we typi- 
cally see that more velocity components are required to fit the absorption. 
Each component has its own physical parameters: central velocity, velocity 
dispersion (rms of thermal and turbulent broadening), ionization, column 
densities and elemental abundances. Prior to its use with quasars, this fit- 
ting method was developed for the ISM, where it represents gas in spatially 
separate clouds. 

Levshakov and co-workers have proposed a different type of model, the 
mesoturbulent model, in which the gas velocities are correlated, and the 
column density per unit velocity is varied to fit the absorption line profiles. 
They assume that the absorption comes from a single region in space, and 
they calculate the distribution of the gas density down the line of sight. To 
simplify the calculations, in early Reverse Monte-Carlo models, they assumed 
that the gas temperature and density were constant along the line of sight, 
which is not appropriate if there are separate discrete clouds of gas with 
differing physical conditions. 

The effects of mesoturbulence on the D/H absorbers towards Q1937-1009 
[0, Q1009+2956 [§7| and Q1718+4807 |Ug] were examined in detail using 



this early model. In the first paper they allowed the Nhi to vary far from 
the observed value (N HI = 7.27 x 10 17 P7|), and consequently they found a 
variety of Nhi, but when the Nhi is held within range, the D/H is 3.3xl0~ 5 , 
exactly the same as with the usual model []62 |. For the second quasar, the 



D/H obtained is again similar to that obtained in the usual way. The results 
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are the same as with the usual model in part because the H and D line 
widths are dominated by thermal and not turbulent motions, and for these 
two quasars the total N HI is not affected, because it is measured from the 
Lyman continuum absorption, which does not depend on velocity. 

Recently they have developed a new model called MCI [f70| , ||101|| appro- 
priate for absorption systems which sample different densities. They now use 
H I and metal ions to solve for two random fields which vary independently 
along the line of sight: the gas density and the peculiar velocities. This 
model allows the temperature, ionization and density to all vary along the 
line of sight. 



The mesoturbulent model of Levshakov et al. |67j and the microturbulent 
Voigt model give the same column densities and other parameters when one 
of the following conditions apply: 1) The line of sight through the absorbing 
gas traverses many correlation lengths. 2) If each velocity in a spectrum 
corresponds to gas at a unique spatial coordinate. 3) The absorbing regions 
are nearly homogeneous, with at most small fluctuations in density or pecu- 
liar velocities, or equivalently, thermal broadening larger than the turbulent 
broadening. 

The Voigt model could give the wrong result when two or more regions 
along the line of sight, with differing physical conditions, give absorption 
at the same velocity. A remarkable and unexpected example of this was 
reported by Kirkman & Tytler ||102|| who found a Lyman limit system which 
comprised five main velocity components. Each component showed both 
C IV and O VI absorption at about the same velocity, but in each of the five 
components, the O VI had a larger velocity dispersion, and hence came from 
different gas than the C IV. While this LLS is much more complex than those 
in which we can see D, this type of velocity structure could be common. 

All authors other than Levshakov and collaborators use standard Voigt 
fitting methods to determine column densities, for several reasons. The Voigt 
method was used, with no well known problems, for many decades to ana- 
lyze absorption in the ISM, and the ISM is well modeled by discrete clouds 
separated in space. The Levshakov et al. [67] methods are more complex. In 
early implementations, Levshokov et al. |)7j made assumptions which are not 
suitable for all absorbers. The current methods require weeks of computer 
time, and in many cases the two methods have given the same results. 

We conclude that, when we have sufficient data, velocity structure is not 
a problem for the absorbers like those now used for D/H. 
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5.5.1 Was the primordial D high but depleted in the absorbers? 

The idea here is that the average BBN D/H was high, and it has been de- 
pleted in the three absorbers which show low D. There are two options: local 
depletion in some regions of the universe, and uniformly global depletion. 
We conclude that there is no known way to deplete D locally, and global 
depletion seems unlikely. 

First we list seven observations which together rule out local depletion, 
including that suggested by Rugers & Hogan ||103|| . 

1. We note that D/H is also low in our Galaxy, and that Galactic chemical 
evolution accounts for the difference from the low primordial D. Hence we 
know of four places where D is low and consistent with a single initial value. 

2. If the BBN D/H was high, let us say ten times larger at 34xl0 -5 , then 
the depletion in all four, widely separated in space, must be by a similar 
factor: Q1937-1009: 0.90 ± 0.02; Q1009+2956: 0.88 ± 0.02; Q0130-4021: 
> 0.80; local ISM in our Galaxy: 0.86 - 0.93, where for the Galaxy alone we 
assume that Galactic chemical evolution reduced the initial D/H by a factor 

of 1.5 3 jgi]. 



3. The quasar absorption systems are large - a few kpc along the line of 
sight |)8|], far larger than can be influenced by a single star or supernovae. 
The gas today in the local ISM is a mixture of gas which was also distributed 
over a similar large volume prior to Galaxy formation. 

4. The abundance of the metals in the quasar cases are very low; too low 
for significant (> 1%) destruction of D in stars |52]. 



5. The quasar absorbers are observed at high redshifts, when the universe 
is too young for low mass stars (< 2 solar masses) to have evolved to a stage 
where they eject copious amounts of gas. 

6. The quasar absorbers are observed at about the time when old stars 
in the halo of our Galaxy were forming. These stars may have formed out of 
gas like that seen in the quasar spectra, but with high density. We expect 
that much of the gas seen in absorption is in the outer halo regions of young 
galaxies, and that some of it was later incorporated into galaxies and halo 
stars. 

7. The ratio of the abundances of Si/C in the quasar absorbers is similar 
to that in old stars in the halo of our Galaxy. This abundance ratio is 
understood as the result of normal chemical evolution. 
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Global destruction of D prior to z = 3, or in the early universe, remains 
a possibility, but it seems contrived. 

Gnedin & Ostriker | |104| | discuss photons from early black holes. Sigl et al. 
[ |105|| show that this mechanism creates 10 times more 3 He than observed, and 
Jedamzik & Fuller [|5^] find the density of gamma ray sources is improbably 
high. 



Holtmann, Kawasaki & Moroi ||106|| , ||107|| showed that particles which 



decay just after BBN might create photons which could photodissociate D. 
With very particular parameters, the other nuclei are not changed, and it 
is possible to get a D/H which is lower than from SBBN with the same f^- 
Hence low D and low Y p can be concordant. An exception is 6 Li which is 
produced with 6 Li/ H ~ 10~ 12 , which is about the level observed in two halo 
stars. There is no conflict with the usual conclusion that most 6 Li is made 
by Galactic cosmic rays prior to star formation, because the observed 6 Li has 
been depleted by an uncertain amount. This scenario has two difficulties: 
Buries (private communication) notes that there would be a conflict with 
the Qb measured in other ways, and it seems unlikely that the hypotheti- 
cal particle has exactly the required parameters to change some abundances 
slightly, within the range of measurement uncertainty, but not catastrophi- 
cally. 

Most conclude that there are no likely ways to destroy or make significant 

D. 



5.5.2 Could the D/H which we observe be too high? 

The answer to this question from Kirshner is, that the D/H could be slightly 
lower than we measure, but not by a large amount. We discuss two possibil- 
ities: measurement problems and biased sampling of the universe. 

First we consider whether the D/H in the quasar absorbers could be less 
than observed. This can readily happen if the D is contaminated by H, 
but a large reduction in D/H is unlikely because the D line widths match 
those expected in Q1937-1009 and Q1009+2956. We do not know how the 
ISM D/H values could be too high, and Galactic chemical evolution requires 
primordial D/H to be larger than that in the ISM, and similar to the low 
value from quasars. Hence it is unlikely that the D/H is much below the 
observed value. 

Second, we consider whether the absorbers seen in the quasar spectra 
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are representative. The absorbers are biased in three ways: they represent 
regions of the universe with well above (100 - 1000 times) the average gas 
density at z — 3, and amongst such high density regions, which are observed 
as Lyman Limit absorption systems, they have relatively low temperatures 
(2 x 10 4 K), and simple quiescent velocity structures. The last two factors are 
necessary to prevent the H absorption from covering up that from D, while 
the high density follows from the high density of neutral H which is needed 
to give detectable neutral D. It is likely that the gas in the absorbers at z = 3 
has by today fallen into a Galaxy, though this is not required because some 
gas will be heated as galaxies form, preventing infall. The low temperatures 
and quiescent velocities argue against violent astrophysical events, and there 
are no reasons to think that the absorbers are any less representative than, 
say, the gas which made up our Galaxy. 

We should also consider whether the quasar absorbers might be unrepre- 
sentative because of inhomogeneous BBN. In this scenario regions with above 
average density will have below average D/H, but the evolution of density 
fluctuations could be such that the low density regions fill more volume ||108| 



47], so that they are more likely to dominate the observed universe today. 
In that scenario the fib derived from the D/H would be below the universal 
average, and the observed (low) value of D/H would be "high" compared to 
expectation for SBBN with the same fib- This scenario will be tested when 
we have observations of many more quasars, 
measurements of D/H towards QSOs. 



5.5.3 Is there spatial variation in D/H towards quasars? 

It seems highly likely that the D is low in the three quasars which show low 
D, and we discussed above why it is hard to imagine how this D could have 
been depleted or created since BBN. Hence we conclude that the low D/H is 
primordial. 

Are there other places where D is high? All quasar spectra are consistent 
with a single low D/H value. The cases which are also consistent with high 
D are readily explained by the expected H contamination. We now explain 
why we have enough data to show that high D must be rare, if it occurs at 
all. 

High D should be much easier to find than low D. Since we have not 
found any examples which are as convincing as those of low D, high D must 
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be very rare. If D were ten times the low value, the D line would be ten times 
stronger for a given Nhi, and could be seen in spectra with ten times lower 
signal to noise, or 100 times fewer photons recorded per A. If such high D/H 
were common, it would have been seen many times in the high resolution, 
but low signal to noise, spectra taken in the 1980's, when the community 
was well aware of the importance of D/H. High D would also have been 
seen frequently in the spectra of about 100 quasars taken with the HIRES 
spectrograph on the Keck telescope. In these spectra, which have relatively 
high signal to noise, high D could be detected in absorption systems which 
have 0.1 of the Nhi needed to detect low D. Such absorbers are about 40 - 
60 times more common than those needed to show low D/H, and hence we 
should have found tens of excellent examples. 



5.5.4 Why is there lingering uncertainty over D? 

Today it is widely agreed that D is low towards a few quasars. There remains 
uncertainty over whether there are also cases of high D, for the following 
reasons: 

• measurements have been made in few places; 

• contamination of D by H looks very similar to D, and resembles high 
D; 

• both the low Y p values reported during the last 25 years, and the 7 Li 
abundance in Spite plateau halo stars, with no correction for depletion, 
imply low fit,, low rj, and high D/H for SBBN; and 

• the first claims were for high D. 

In most cases, the apparent conflicts over D/H values concern whether 
the absorption near the expected position of D is mostly D or mostly H. 
Steigman ||109|| and all observational papers discussed this contamination of 
D by H. 

Carswell et al. pO] noted that contamination was likely in Q0014+813 



and hence the D/H could be well below the upper limit. Songalia et al. |7T 
stated: "because in any single instance we can not rule out the possibility of 
a chance H contamination at exactly the D offset, this result [the high D/H] 
should be considered as an upper limit until further observations of other 
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systems are made." Buries et al. [|72j showed that Q0014+813 is strongly 
contaminated, does not give a useful D/H limit. For Q1718+4807 we |)8 



and Levshokov, Kegel & Takahara | |100| | have argued that contamination is 
again likely. 

There are many reasons why contamination is extremely common: 

• H absorption looks just like that from D, 

• H is 30,000 times more common, 

• spectra of about 50 quasars are needed to find one example of relatively 
uncontaminated D, 

• high signal to noise spectra are needed to determine if we are seeing H 
or D, and 

• these spectra should cover all of the Lyman series and metal lines, 
because we need all possible information. 

When H contaminates D, the resulting D/H will be too high. 

It is essential to distinguish between upper limits and measurements. 
There are only two measurements (Q1937-1009 and Q1009+2956). They are 
measurements because we were able to show that the D absorption line has 
the expected width for D. All other cases are upper limits, and there is no 
observational reason why the D/H should be at the value of the limit. In 
many cases, all of the D can be H, and hence and D/H = is an equally 
good conclusion from the data. 

Only about 2% of QSOs at z ~ 3 have one absorption systems simple 
enough to show D. All the rest give no useful information on D/H. Typically, 
they do not have enough H to show D, or there is no flux left at the position 
of D. In such cases the spectra are consistent with high, or very high, D/H, 
but it is incorrect to conclude that D/H could be high in ~ 98% of abosrption 
systems because these systems are not suitable to rule out high D/H. Rather, 
we should concentrate on the few systems which could rule out both high 
and low D/H. 

We will continue to find cases like Q1718+4807 which are consistent with 
both low and high D/H. As we examine more QSOs we will find some cases of 
contamination which look exactly like D, even in the best spectra, by chance. 
But by that time we will have enough data to understand the statistics of 
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contamination. We will know the distribution function of the contaminat- 
ing columns and velocities, which we do not know today because the D/H 
absorbers are a rare and special subset of all Lyman limit absorbers. When 
absorbers are contaminated we will find a different D /H in each case, because 
the Nhi, velocity and width of the contaminating H are random variables. 
But we will be able to predict the frequency of seeing each type of contami- 
nation. If there is a single primordial D /H then we should find many quasars 
which all show this value, with a tail of others showing apparently more D/H, 
because of contamination. We will be able to predict this tail, or alterna- 
tively, to correct individual D/H for the likely level of contamination. When 
we attempted to correct for contamination in the past f53|, ||110|| , f68fl , we 
used the statistics of H I in the Lya forest because we do not have equivalent 
data about the H I near to the special LLS which are simple enough to show 
D/H. Such data will accumulate at about the same rate as do measurements 
of D/H, since we can look for fake D which is shifted to the red (not blue) 
side of the H I. 

There are large differences in the reliability and credibility of different 
claimed measurements of D/H in quasar spectra, and hence much is missed 
if all measurements are treated equally. It also takes time for the commu- 



nity to criticize and absorb the new results. Early claims of high D/H ||103 



111|| in Q00 14+81 18 are still cited in a few recent papers, after later mea- 



surements |72|] with better data, have shown that this quasar gives no useful 
information, and that the high D/H came from a "spike" in the data which 
was unfortunately an artifact of the data reduction. 

In summary, the lack of high quality spectra, which complicates assess- 
ment of contamination by H, is the main reasons why there remains uncer- 
tainty over whether some absorbers contain high D. 



5.5.5 Why we believe that the D/H is Primordial 

Here we review why we believe that the low D/H is primordial. These ar- 
guments are best made without reference to the other nuclei made in BBN, 
because we wish to use the abundances of these nuclei to test SBBN theory. 

• D/H is known to be low in four widely separated locations: towards 
three quasars, and in the ISM of our Galaxy. 
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The extraction of D/H from quasar spectra is extremely direct, except 
for corrections for contamination by H, which make D/H look too large. 

Since contamination is common, all data are consistent with low D/H, 
and no data require high D/H. 

High D/H is rare, or non-existent, because it should be easy to see in 
many existing spectra, but we have no secure examples. 

The low D/H in the quasars, pre-solar system and in the ISM today 
are all consistent with Galactic chemical evolution. 

The quasar absorption systems are large - many kpc across, as was the 
initial volume of gas which collapsed to make our Galaxy. 

The abundance of the metals in the quasar cases are very low, and 
much too low for significant (> 1%) destruction of D in stars. 

The quasar absorbers are observed at high redshifts, when the universe 
is too young for low mass stars to have evolved to a stage where they 
eject copious amounts of gas. 

The ratio of the abundances of Si/C in the absorbers is normal for 
old stars in the halo of our galaxy, indicating that these elements were 
made in normal stars. 

In the quasar absorbers, the temperatures and velocities are low, which 
argues against violent events immediately prior to the absorption. 

If BBN D/H were high, the hypothetical destruction of D would have 
to reduce D/H by similar large amounts in all four places. 

The above observations make local destruction of D unlikely. 

There are no known processes which can make or destroy significant D. 

Global destruction of D by photodissociation in the early universe re- 
quires very specific properties for a hypothetical particle, and is limited 
by other measures of 
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5.5.6 Conclusions from D/H from quasars 



Most agree that D is providing the most accurate 77 value (JTTJ , although some 
have one remaining objection, that there might also be quasar absorbers 
which show high values of D/H ||112j| , J7J. 

The D/H from our group (Buries & Tytler |)7|], Q, [f^Q, together with 



over 50 years of theoretical work and laboratory measurements of reaction 
rates, leads to the following values for cosmological parameters (unlike most 
errors quoted in this review, which are the usual la values, the following are 
quoted with 95% confidence intervals): 

• D/H = 3.4 ± 0.5 x 10~ 5 (measured in quasar spectra) 

• 77 = 5.1 d= 0.5 x 10~ 10 (from BBN and D/H) 

• Y p = 0.246 ± 0.0014 (from BBN and D/H) 

• 7 Li/H = 3.5lJi x 10" 10 (from BBN and D/H) 

• 411 photons cm -3 (from the CMB temperature) 

• p b = 3.6 ± 0.4 x 10" 31 gcm- 3 (from CMB and 77 ) 

• tt b h 2 = 0.019 ± 0.0024 (from the critical density p c ) 

• N v < 3.20 (from BBN, D/H and Y p data). 

If we accept that D/H is the most accurate measure of 77, then observa- 
tions of the other elements have two main roles. First, they show that the 
BBN framework is approximately correct. Second, the differences between 
the observed and predicted primordial abundances teach us about subsequent 
astrophysical processes. Recent measurements of 4 He agree with the pre- 
dictions. It appears that some 7 Li has been destroyed in halo stars 1 1 1 3 1 , 
and 3 He is both created and destroyed in stars. 



6 Helium 

The high abundance of 4 He allows accurate measurements in many locations. 
However, 4 He is also produced by stars, and since such high accuracy is 
required, the primordial abundance is best measured in locations with the 
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least amounts of stellar production. High accuracy is desired, since D/H 
predicts Y p to within 0.0014 (SY P /Y P = 0.006, 95% confidence), which is well 
beyond the typical accuracy of astronomical abundance determinations. In 
the local ISM, the amount of 4 He from stars is about Y = 0.01 — 0.04; much 
less than Y p , but ten times the desired accuracy for Y p . 

Helium has been seen in the intergalactic medium, where Carbon abun- 
dances are < 0.01 solar, and possibly zero in much of the volume. Strong 
absorption is seen from the He II Lya line at 304A in the redshifted spectra 
of quasars [|114|1 , however it is difficult to obtain an abundance from these 



measurements, because nearly all He is He III which is unobservable, and we 
do not know the ratio He II/He to within an order of magnitude. However, 
the strength of the He II absorption does mean that there is abundant He in 
the intergalactic gas [ 115| |, which has very low metal abundances, which is 
consistent with BBN, and probably not with a stellar origin for the 4 He. 

The best estimates of the primordial abundance of He are from ionized gas 
surrounding hot young stars (H II regions) in small galaxies. The two galaxies 
with the lowest abundances have 1/55 and 1/43 of the solar abundance. The 
4 He and H abundances come from the strengths of the emission lines which 
are excited by photons from near by hot stars. 

Values for Y p from these extragalactic H II regions have been reported 
with small errors for more than 25 years, e.g.: 



Y p = 0.216 ± 0.02 [34 



Y p = 0.230 ± 0.004 gT6 



Y p = 0.234 ± 0.008 [|TT 



Y p = 0.236 ± 0.005 [|TT 



• Y p = 0.228 ± 0.005 [|TT 

• Y p = 0.234 ± 0.002 ± 0.005 |6f . (random, and systematic errors) 

• Y p = 0.246 ± 0.0014 (95% prediction from low D/H and SBBN). 

These values are lower than the value now predicted by low quasar D/H 
and they appear incompatible, because of the small errors. However Skillman 
et al. [ |120|| argued that errors could be much larger than quoted, allowing 
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Y p < 0.252, and Pagel ||121|| (and private communication 1994) agreed this 
was possible. 

The measurement of Y p involves three steps. Emission line flux ratios 
must be measured to high accuracy, which requires good detector linearity 
and flux calibration, and corrections for reddening and stellar He I absorp- 
tion. These fluxes must be converted to an abundance, which requires correc- 
tion for collisional ionization and neutral He. Correction for unseen neutral 
He depends on the spectral senegy distribution adopted for the ionizing radi- 
ation and might change Y p by 1 — 2 percent. Then the primordial abundance 
must be deduced from the Y values in different galaxies. 

Izotov, Thuan & Lipovetsky ||122 



12311 have been pursuing a major 



observational program to improve the determination of Y p . They have found 
many more low metallicity galaxies and have been reporting consistently 
higher Y p values, most recently in their clear and persuasive paper [T|]: 

• Y p = 0.244 ± 0.002 from regression with O/H and 

• Y p = 0.245 ± 0.001 from regression with N/H. 

The four main reasons why these values are higher are as follows, in order 
|125|| , (Skillman, and Thuan personal communication 



124 



of importance | 
1998). 

1. When stellar He I absorption lines underlying He emission lines are 



not recognized, the derived Y p is too low. This is a important for IZwl8 [|124 



which has the lowest metallicity and hence great weight in the derivation of 
Y p , and perhaps for many other galaxies. 

2. The emission line fluxes must be corrected for collisional excitation 
from the metastable level. At low abundances, which correlate with high 
temperatures, these corrections can be several percent. The amount of cor- 
rection depends on the density. There are no robust ways to measure these 
densities, and differing methods, used by different groups, give systematically 
different results. Izotov and Thuan ||124|| solve for the He II density, while 
Olive, Skillman and Steigman [PJ use an electron density from the S II lines. 



3. Izotov & Thuan ||124|| have spectra which show weaker lines, and they 
use the five brightest He lines, while Olive et al. 
6678. 

4. Izotov & Thuan [fT 



46| usually use only Hel 



1241 correct for fluorescent enhancement, which 



increases the Y values from for a few galaxies. 
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For these reasons Izotov & Thuan [ |124|| obtain higher Y values for indi- 
vidual galaxies which have also been observed by Olive, Skillman & Steigman 
PJ, and Izotov & Thuan \ VZA \ find a shallower slope for the regression to 
zero metal abundance (see |ll| Fig 6). And most importantly, using higher 
quality Keck telescope spectra, they obtain high Y p = 0.2452±0.0015 (random 



errors), from the two galaxies with the lowest metal abundances ||126| 



These measurement difficulties, combined with the recent improvements, 
lead most to conclude that the Y p is in accord with the SBBN. The Izotov & 
Thuan fl| values are very close to the low D/H predictions, while the lower 
Y p quoted by Olive |7j, 0.238 ± 0.002 ± 0.005, is also consistent when the 
systematic error is used. 

It is clear that the systematic errors associated with the Y p estimates have 
often been underestimated in the past, and we propose that this is still the 
case, since two methods of analyzing the same Helium line fluxes give results 
which differ by more than the quoted systematic errors. While the Izotov & 
Thuan |124fl method has advantages, we do not know why the method used 
by Olive, Skillman & Steigman j|6| should give incorrect answers. Hence the 
systematic error should be larger than the differences in the results: 0.007 
using the most recent values, or 0.011 using earlier results. 



7 3 He 

The primordial abundance of 3 He has not been measured. This is most 
unfortunate, since it is nearly as sensitive as D to the baryon density during 
BBN. 3 He is harder to measure than D because the difference in wavelength 
of 3 He and 4 He lines is smaller than for D, and the Lyman series lines of 
He II, main absorption lines of He in the IGM, are in the far ultraviolet at 
228 - 304A which is hard to observe because of absorption in the Lyman 
continuum of H I at < 9 12 A. 

Rood, Steigman & Tinsley [|127]| argued that it was unlikely that 3 He could 



be used to supplement cosmological information from D because low mass 
stars should make a lot of 3 He, increasing the current ISM value to well above 
that in the pre-solar system ISM, and in potential conflict with observations 
at that time. This conflict has been confirmed. Measurements do show 
enhanced 3 He in Planetary nebulae, as expected from the production in the 
associated low mass stars, but this is not reflected in the ISM as a whole. 
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The pre-solar and current 3 He abundances are similar [|128 
with expectation ||129|| , [|130|| , for unknown reasons. 
It was suggested ( [fL31|| , see review by Hata et al. 



in contradiction 



13211) that the uncer- 



tainty over the amount of destruction of D could be circumvented using the 
sum of the abundances of D + 3 He, since the destroyed D should become 3 He, 
and 3 He is relatively hard to destroy. The primordial D + 3 He should then 
be < the same sum observed today, as more 3 He is made in stars over time. 
However, there are two problems with this scenario. First, the 3 He should 
increase over time, which it does not, implying that some stars destroy 3 He, 
and second, the 3 He abundance should be about constant in the ISM today, 
which it appeared not to be in early data 
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Hence, just prior to the 
measurement of D in quasars, most concluded that D + 3 He in the Galaxy 
does not provide secure cosmological information [ |132j| , and summaries by 

ra, irron. 



Balser et al. 



|134|| report on a 14 year program to measure 3 He in the 
Galactic H II regions. Using models for the gas density structure, they find 
an average 3 He/H = 1.6±0.5 x 10~ 5 for a sub-sample of seven simple nebulae. 
No variation is seen with Oxygen abundance over a factor of ten, and there 

This value may represent the average in the ISM today, 



is little scatter 135 



but it is not known how to use this to obtain primordial abundances. 

These measurements are relevant to stellar nucleosynthesis and Galactic 
chemical evolution, and are consistent with a cosmological origin for the 3 He, 
but we suggest that gas with much lower metal abundances will need to be 
observed to derive a secure primordial abundance for 3 He. 



8 Lithium 



Lithium is observed in the solar system, the atmospheres of a wide variety 
of stars and in the ISM. Arnould & Forestini 



l i L36| ] review light nuclei abun- 
dances in a variety of stars and related stellar and interstellar processes, while 
halo stars are reviewed by pS| , |137|| , [ |113| | and |4T5|. 

Old halo stars which formed from gas which had low iron abundances 
show approximately constant 7 Li/H ~ 1.6 x 10~ 10 and little variation with 
iron abundance or surface temperature from 5600 - 6300 K. The lack of 
variation amongst these "Spite plateau" stars ||138|| (references in ||J) shows 
that their 7 Li is close to primordial. 



31 



Since the halo stars formed about ten times more 7 Li has been produced 
in the inner Galaxy. Abundances of 7 Li/H ~ 10~ 9 are common, although 
some stars show more, presumably because they make 7 Li. Stars typically 
destroy 7 Li when they evolve, accounting for the low abundances, < 10 -11 , 
in evolved stars. Stars with deeper convection zones, such as halo stars with 
lower surface temperatures, show less 7 Li, because they have burnt it in their 
interiors. 

Here, and in the next section on 6 Li, we will the following topics: 

• measurement of current surface abundances on the Spite plateau, 

• change in 7 Li with iron abundance, 

• creation of 7 Li and 6 Li after BBN and prior to halo star formation, 

• depletion of these nuclei in the halo stars, 

• stars with differing 7 Li, and 

• gravitational settling. 



The recent homogeneous data on 22 halo stars with a narrow range of 
temperature on the "Spite plateau" have very small random errors and show 
that most (not all) stars with similar surface properties have the same 7 Li/H 



113]. Earlier data showed more scatter, which some considered real (refer- 
ences in [|139|| ), and hence evidence of depletion. 

The Ryan, Norris, & Beers [|113|| sample shows a clear increase of 7 Li 
with iron abundance, as had been found earlier. This trend appears to be 
real, because the data and stellar atmosphere models used to derive the 7 Li 
abundance do not depend on metallicity. But it was not found by Bonifa- 
cio & Molaro ] [140|| , perhaps because of larger scatter in temperatures and 
iron abundance. This trend is not understood, and there are several possible 
explanations. It may have been established in the gas from which the stars 
formed, perhaps from cosmic rays in the ISM, or from AGB stars. Alter- 
natively, we speculate that it might instead relate to depletion of the 7 Li in 
the stars. In either case, the BBN 7 Li will be different from that observed: 
smaller if the 7 Li was created prior to the star formation, and higher if the 
trend is connected to destruction in the stars. More on this below. 
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Creation of 7 Li in the ISM by cosmic ray spallation prior to the formation 
of the halo stars is limited to 10 - 20% because Be would also be enhanced 
by this process ||113|| , ]7). 

A clear summary of arguments for and against significant depletion is 
given by Cayrel [|8|. There are two main reasons why depletion is believed to 



be small: the negligible dispersion in 7 Li for most halo stars on the plateau, 
and the presence of 6 Li. The main arguments for depletion are that it is 
expected, it clearly occurs in some stars, some halo stars on the plateau show 
differing abundances, and star in the globular cluster M92 which should have 
similar ages, composition and structure, show a factor of two range in 7 Li. 

Different depletion mechanisms include mixing induced by rotation or 
gravity waves, mass loss in stellar winds and gravitational settling. Some 
models predict either variation from star to star, or trends with temperature, 
which are not seen for the stars on the plateau. For example, the rotationally 
induced mixing model implies that stars with different angular momentum 
histories will today show different 7 Li. Ryan, Norris & Beers [|113|1 find that 
the small scatter in their data, especially after the removal of the correlation 
with the iron abundance, limits the mean depletion in these models to < 30%, 
much less than the factor of two needed to make 7 Li agree exactly with the 
predicted abundance from low D/H. 

Some stars which should lie on the plateau have very low 7 Li, while others 
show a range of abundances (see ref. in ||113|| ). Differences are also seen 



between halo field stars [ |113|1 and stars in the globular cluster M92 [|141 



143 ], which show a factor of two spread in 7 Li. These observations are not 
understood. 

Gravitational settling (diffusion) of heavier elements reduces the 7 Li in 
the atmospheres of stars. However, the depletion should be most in the 
hottest (highest mass) stars, which is not seen, and not understood. Vauclair 
& Charbonnel [144] proposed that small stellar winds might be balancing 
the settling. Vauclair & Charbonnel | |145| noted that the peak abundances 
inside the stars are independent of both mass and iron abundance. Normal 
stellar models predict that these peak abundances will not be seen in the 
stellar atmospheres, because convection does not reach this far down into 
the stars. However they point out that if some mechanism does mix gas 
from the 7 Li peak zone into the bottom of the convection zone then the 
stars on the plateau would have similar abundances as observed. Assuming 
that the observed abundances are those from the peaks inside the stars, they 
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find that the initial abundance in the stars was 7 Li/H = 2.2 ± 0.6 x 10~ 10 , 
without free parameters, which is still below but statistically consistent with 
the prediction from low D of 3.5to'g xl0~ 10 . 



8.1 Primordial 7 Li 



Ryan, Norris and Beers [|113|1 conclude 7 Li/H ~ 10~ 10 , with small random 
errors and three sources of systematic error, each up to a factor of 1.3, from 
the effective temperatures, stellar atmospheres and enhancement prior to 
star formation. Bonifacio & Molaro |I30| found 7 Li/H = 1.73 ± 0.05 ± 0.2 



xlO 10 . These abundances are both below the value of 3.5l£;£ x lO" 10 (95%) 
from BBN and our D/H, but unlike |113|, we feel they are not inconsistent 



given the quoted systematic errors, the lack of understanding of depletion, 
and the variation amongst similar stars. We do not know how to estimate 
the systematic errors connected with these issues. Given the comparative 
simplicity of D/H, we prefer to use it and SBBN, and we stick with our 
earlier suggestion ]53| that 7 Li in the Spite plateau halo stars is depleted 
by about a factor of two. Most, but not all agree that this is reasonable. 
Depletion by much larger factors, which was discussed a few years back, 
is now our of favor because of improved models. Improved modelling of 
rotational mixing, has lead to better fits to high metal abundance (population 
I) stars, which can be applied to halo (population II) stars, while the initial 
rotation rates of the halo stars may be lower than was assumed (Deliyannis 
private communication) . 

In summary, both the data and theory tells us that the 7 Li on the Spite 
plateau is not exactly the primordial value. The correction is probably small, 
less than a factor of two, but we do not yet know its value. 

If we are to attain a primordial 7 Li abundance we must either (1) un- 
derstand why its abundance varies from star to star, and learn to make 
quantitative predictions of the level of depletion, or (2) make measurements 
in relatively unprocessed gas. 

We are optimistic that primordial 7 Li will be measured to high precision. 
Compared to D and He, the observations are simple: 15 - 20 mA absorption 
lines in relatively empty spectra of often bright stars (V=ll). The best 
data have small errors. We anticipate that further studies will determine 
the amount of 7 Li produced prior to the formation of the stars, and the 
subsequent depletion in these stars. The possible increase in 7 Li with iron 
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abundance is a clue, as are the 6 Li, Be and B abundances in the same stars. 



8.2 6 Li 



The primordial 6 Li abundance has not been observed, but 6 Li/H has been 
measured in two stars on the Spite plateau. The abundance is well below that 
expected from SBBN, but 7 Li is used to help determine the primordial BBN 
6 Li abundance in two ways. First, the presence of 6 Li limits the amount 
of destruction of 7 Li, because 6 Li is more fragile than 7 Li. Second, if the 
observed 6 Li was made prior to the formation of the stars, then some, perhaps 
much 7 Li [ |139|| , may have been made by similar processes. The first point 



is often presented as evidence that the 7 Li on the Spite plateau is close to 
primordial (e.g. less that a factor of two depletion, according to [f78|), but 
the second point is cause for caution. 

6 Li has been detected in only two stars on the Spite plateau, because the 
absorption line at 6707.97A is weak and fully blended with 7 Li at 6707. 8lA. 
This is a difficult observation. The 6 Li makes the absorption line slightly 
asymmetric, and this is detected using models of the line broadening, which 
are tested on other absorption lines which are expected to have similar profiles 
because they arise in the same layers of the stellar atmosphere. Following the 
impressive first detection by [ |146|| and [147|, and Cayrel et al. ||148|| report 
6 Li/ 7 Li= 0.052 ± 0.019 in HD84937, while Smith, Lambert & Nissen 
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report 6 Li/ 7 Li= 0.06 ± 0.03 in BD+26 3578. It is not known whether these 
detections are representative of halo stars on the Spite plateau. Most assume 
that they are, but they could be above normal, perhaps by a lot; Smith et 
al. report 6 Li/ 7 Li= 0.00 ± 0.03 for six other stars. 
The SBBN makes 6 Li/H 



10 



-13.9 



1491 , 1150 



using the rj from D/H, 
which is 500 times less than the measured abundance of 7 x 10~ 12 in the two 
halo stars. The SBBN isotope ratio is 6 Li/ 7 Li= 3 xlO -5 , a factor of 2000 less 
than observed in these two stars. This is not considered a contradiction with 
SBBN, because 6 Li, and some 7 Li at the same time, can be made elsewhere. 

The 6 Li is usually assumed to have been present in the gas when the stars 
formed, but it could be made later, e.g. when cosmic rays strike the star or 
in stellar flares [ |141| |. Production by cosmic rays in the ISM prior to the star 
formation is most favored | 142 |. With this assumption, the effects on 7 Li 
can be calculated in two steps. First, determine the ratio of 6 Li/ 7 Li in the 
production process (the production ratio). Second, correct for the depletion 
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of 6 Li in the stars to determine the initial abundance of 6 Li. The amount of 
7 Li produced along with the initial 6 Li is then specified. 

Cosmic rays in the early ISM could have made 6 Li and some 7 Li prior to 
the formation of the Spite plateau halo stars. The production ratio depends 
on the reaction and energies (e.g. ||139|| ). Two reactions of cosmic rays in 
the ISM are considered to produce 6 Li. Smith, Lambert & Nissen [|139|| find 
that 6 Li/Be ratios imply that most 6 Li was made in a — a fusion reactions, 
rather than in spallation (e.g. O + p — > 6 Li) which is favored by ||142|| and 
150]. The production ratio is 6 Li/ - 7 Li ~ 2 for the a — a reaction. 



Standard stellar models ||141|| predict that much of the initial 6 Li will 
have been destroyed in the stars. The more that was destroyed, the more 
6 Li and non-BBN 7 Li should have been in the initial gas to give the observed 
abundances. Depending on the destruction mechanism, the destruction of 
6 Li may also destroy 7 Li, but this is usually ignored. 

When we choose the amount of depletion of 6 Li, we fix the amount present 
when the stars formed. If the 6 Li has been depleted by a large factor, ~ 100, 
then the stars would have begun with 6 Li/ 7 Li similar to the production ratio, 
and essentially all of the 7 Li would be non-primordial ||139|| , which is an 
unusual conclusion! 

Ryan, Norris and Beers ||113|| assume that 50% of the 6 Li and none of the 
7 Li was destroyed, and use a production ratio of 1.5 to conclude that the 
BBN 7 Li was 0.84 of that now in the stars. Since nearly all observations of 
Li are made at low resolution, the 6 Li and 7 Li lines are not resolved, they 
correct for the 6 Li . If the two stars with observed 6 Li are normal, then the 
BBN 7 Li is about 79% of the observed Li absorption. 

Many other papers discuss this topic. Olive & Fields |14[ give a summary. 
Cayrel et al. ||148|| use models for the formation of Li, Be and B and calculate 



the expected abundance of 6 Li when the star formed, and find that the 
observed abundance implies little depletion of 6 Li, and a 7 Li depletion of 
less than 25%. Vangiono-Flam et al. | |150[ | also argue that 6 Li is not much 
depleted, and find that its BBN abundance, extrapolated back to before the 
production by spallation, is compatible with a BBN abundance of 3 x 10~ 13 
- 5.6 x 10~ 14 . 

All eagerly await the measurement of 6 Li, together with Beryllium, in 
more stars. 
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9 Beryllium 



The primordial abundance of Beryllium has not been observed, 
duction in SBBN is 9 Be/H < 1(T 17 |U9 
the observed level. Inhomogeneous BBN allows much higher abundances, 



The pro- 
150|1 , orders of magnitude below 



possibly approaching detection ||149 



Be is observed. It is created in the ISM when cosmic rays strike C, N and 
O nuclei, and it is destroyed in stars. It is difficult to use Be to constrain 
the cosmic ray production of Li because the production ratio is highly model 
dependent [ |151|| . Beryllium is observed in the atmospheres of halo stars, 
including those on the Spite plateau. Boesgaard et al. [ |152| [ have found that 
Be increases with Iron, and that Be increases 8 times faster than Oxygen, a 
rate consistent with cosmic ray creation. There is some evidence for a spread 
in Be as a given Fe/H, but no sign of a primordial plateau, down to Be/H 

= 10 -13.5_ ' 



10 Are the different nuclei concordant or is 
there a crisis? 



Nearly everyone believes that the primordial abundances are consistent with 

171 



BBN (e.g. ITU, 0, § 



T54f 



questions about the measurements. 



but there are many lingering 
The reader will readily detect the two 
attitudes described by Audouze |5H}[ : "optimistic" , and "agnostic and perhaps 
heretical" in many papers. Each of us tends to adopt differing attitudes for 
each nucleus and astrophysical processes. This review favors D/H because it 
is simple and familiar. 
Steigman 



noted that there was "a hint of an emerging crisis" be- 
cause the 4 He abundances appeared to be lower than expected using the rj 
from the other nuclei, but he recommended much more careful study of the 
uncertainty in BBN predictions, chemical evolution, and observational un- 
certainties including systematic effects. Hata et al. [ 155 and Steigman [156 
stated that "there is a conflict" , referring to the differences in r\ implied by 
low D and low Y p values. 

Whether or not there is a crisis depends on the confidence assigned to 
the answers to three questions: 
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• Is primordial D /H low everywhere, or are there also a some places with 
high values? 

• Is Y p low, high, or uncertain? 

• Has 7 Li in halo stars been depleted by a factor of two? 

Some combinations of answers are not consistent with SBBN. Recent 
data make low D/H seem secure in three quasars plus the ISM, hence the 
issue is whether there are also other places with high primordial D. Low 
D/H is compatible with high Y p and depleted 7 Li, but not with low Y p or 
undepleted 7 Li. High D is compatible with low Y p and undepleted 7 Li, but it 
is incompatible with the three sites which show low D/H and with Galactic 
chemical evolution. A factor of ten D depletion would be required in all 
four places. Low Y p is compatible with undepleted 7 Li and high D, but is 
incompatible with the low D. 

A good case has been made for high Y p , explanations have been given 
why earlier results gave lower values, and the uncertainty appears to be 
larger than quoted. Hence D and Y p are in agreement. 

The 7 Li observed in stars on the Spite plateau is lower than values consis- 
tent with low D. Depletion might provide an explanation, but the amount of 
depletion and the dominant mechanism are not known. The lack of scatter 



implies little depletion, less than expected, which [113], [4£] conclude is not 



sufficient to match low D. Bonifacio & Molaro ||140|| find a higher 7 Li, but 



still below the level required to match low D without depletion. 



11 Non-standard BBN 



The many different forms of non-standard BBN have been reviewed by Coles 
& Lucchin |22[] and Jedamzik ||157|| . Much work has been devoted to inhomo- 



geneous baryon distributions during BBN, additional relativistic particles, 
decaying particles, large neutrino chemical potentials (e.g. | |158[ |), sterile 
neutrinos (e.g. ||159|| ), magnetic fields (e.g. ||160|| ), anti- matter domains (e.g. 
[|161|| ), and alternative theories of gravity (e.g. 
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11.1 Inhomogeneous BBN 

Following early discussion of inhomogeneous BBN (IBBN) by Epstein & Lat- 
timer ||163|| and Hogan ||164| , many detailed studies of different types of inho- 



moegeneity have been published. Malaney & Mathews |]165[| and Kainulainen, 
Kurki-Suonio, & Sihvola ||154j| give reviews. 



IBBN has been discussed to allow larger than standard BBN, to allow 
differing values of D/H in the universe, and to reconcile low Y p with low D/H 
values. 

One exciting goal of this work was to determine whether inhomogeneity 
could give the observed abundances with Q b much larger than the usual value, 
and perhaps large enough to account for all gravitating matter, without the 
need for non-baryonic dark matter (e.g. ||164| , ||104|| , [|108|| ). The best upper 



limit on Of, comes from the lowest observed D/H, which until recently was in 
the ISM. In standard BBN, a higher is ruled out because BBN would make 
less than the observed ISM D/H, and no other way to make D is known. In 
IBBN the D/H in the ISM comes from low density regions, allowing a higher 
average density. The current observations, with some exceptions, fit SBBN 
well, and hence IBBN allows only a slight increase in Q b . 

Inhomogenieties can be imagined over a wide range of distance scales. 
The smallest scales, < 10 -5 pc, mix prior to BBN, leaving homogeneous 
SBBN. Small scales mix during BBN. Intermediate scales which mix after 
BBN give abundances which are constant in space today, but the abun- 
dances are different from SBBN with the same fib- Extra D would be made 
in regions with low density during BBN, giving enhanced D/H everywhere 
today. Large scales (> 1 kpc) may have avoided mixing, and could give dif- 
ferent D/H in different locations today. The near isotropy of the CMB limits 
inhomogenieties to < 1 Mpc. 

Jedamzik & Fuller |^7[ found it difficult to match observed abundances 
of 7 Li with large scale primordial isocurvature baryon number fluctuations. 
Most overly dense regions of the universe with masses greater than the local 
baryon Jeans mass would have to collapse (to prevent observation of the 
7 Li which is overproduced) and smaller scale fluctuations would have to be 
absent or suppressed. Gnedin, Ostriker & Rees [ [L66|| and Copi, Olive & 



Schramm ||167|| reached similar conclusions. Copi, Olive & Schramm ||168 



also showed that large scale (>>lMpc) isocurvature perturbations conflict 
with the smoothness of the CMB, but do not rule out inhomogeneity [0. 
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Kainulainen, Kurki-Suonio, & Sihvola review IBBN. The fib can be 
higher than in SBBN provided the distance scale of the baryon inhomogeneity 
is near to optimal to maximize neutron diffusion effects. The distance scale 
expected for inhomogeneities arising in the electroweak transition are too 
small (10~ 6 to 10~ 3 pc today) to have major effects, although not below the 
accuracy of BBN abundance calculations. QCD inhomogeneities are not so 
limited. However, a low D/H < 5 xIO -5 still requires Y p > 0.240 even in 
IBBN, which helps reconcile low D/H and low Y p measurements, especially 
when we accept that the errors on Y p are larger than quoted. 

Rehm & Jedamzik ||f 6I|| studied BBN in the presence of anti-matter do- 



mains. Annihilation is preferentially on neutrons, and in a limiting case the 
resulting universe is without light nuclei, in violation of the measured abun- 
dances. With small amount of anti- matter, both the low Y p and low D/H 
measurements are matched. 

Early IBBN results looked promising. Today it appears that the scales are 
too small to have major effects, and measurements of primordial abundances, 
especially upper limits on 7 Li, with modest depletion (< factor of two), are 
usually used to give limits on the inhomogeneity, rather than to argue that 
inhomogeneity helps explain discordant data or allows different conclusions 
about fib- 

11.2 The number of Relativistic Particles and their 
Decays 

The main idea here is that the 4 He abundance depends on the number of 
relativistic particles during BBN. Extra particles, such as neutrinos or su- 
persymmetric particles, which are relativistic during BBN, lead to faster 
expansion, larger n/p and a larger Y p . 



Steigman, Schramm & Gunn ||I69|| calculated that BBN limited the num- 



ber of families to N u < 5 to match the 4 He abundance. The range allowed 
by SBBN and laboratory measurements have both narrowed over the years 
and agree well today HTTPl Q. A recent update f7| gives N v < 3.20 (95%) 



from SBBN, although a larger range is obtained if a wider variety of mea- 
sured abundances are accepted [|I7I|| . March- Russell et al. ||158|| note that 



additional relativistic degrees of freedom are allowed if there is a large com- 
pensating asymmetry in the electron neutrino number. 
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Shi, Fuller & Abazajian ||159|| follow the time evolution of a lepton num- 
ber asymmetry arising from active - sterile neutrino transformations during 
BBN. For v e mixing with u s , Y p was allowed to change from -1% to +9%, 
while v T or mixing with v s allowed -2% to +5%. Hence the Y p predicted 
by low D/H in SBBN could be lowered to 0.241, which is between the high 
and low measurements. 

For many years past, observations suggested that Y p was smaller than 
expected for the low D/H in the ISM and now QSOs. It is hard to make Y p 
lower, since this requires fewer, not more, particles than in SBBN. Holtmann 
| p. 071 proposed decays of neutrinos, but this is nearly ruled out 
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by the Kamiokande results on atmospheric neutrinos ||154| . 

Lindley |fL72|| found that massive particles decaying into photons must 
have lifetimes in excess of a few thousand seconds, to avoid the destruction of 
BBN D. Audouze, Lindley & Silk [1173] noted that such radiative decays could 
photodisintegrate 4 He and make D and 3 He, removing the upper bound on 
Qb- Other references are given by Holtmann et al. [|107fl , who discuss weakly 
interacting massive (100 GeV) particles which decay of order 10 6 s after 
BBN. The authors give limits on the abundance and lifetimes of gravitinos 
and neutralinos, for a wide range of light nuclei primordial abundances. 

Kohri & Yokoyama ||174|| give limits on the mass fraction in primordial 
black holes with masses 10 8 — 3 x 10 10 g which evaporate during BBN and 
change the abundances. 

Lopez-Suarez & Canal |175| combine inhomogeneous nucleosynthesis and 
particles which decay at a late time to reassess the limits on Qb- They find 



parameters which allow f2fe< 0.13 



0.18/ifo 2 



(h 70 is the Hubble constant in 



units of 70 km s Mpc ). Such high fib might appear to remove the need 
for non-baryonic dark matter, but there would then be conflicts with other 
measures of fib, especially the baryon fraction in clusters of galaxies, if all 
those baryons were observable today. 



12 Cosmo logical Baryon Density 

The measurement of the baryon density is now a highly active area of re- 
search. In the coming years, we anticipate that higher accuracy measure- 
ments of the baryon density, from the CMB, clusters of galaxies, and the 



Lya forest, will give a new rigorous test of BBN |TTJ. This test can be 
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viewed from two directions. First, we can use the baryon density to fix the 
last free parameter in BBN, and second, we can compare the different baryon 
density measurements, which should be identical if SBBN is correct, and all 
baryons are counted in the measurements made at later times. 

In addition to BBN, the baryon density is measured in four ways: in 
the IGM, in clusters of galaxies, using simulations of galaxy formation, and 
directly from the CMB. All agree with the value from SBBN using low D/H, 
but today they are each about an order of magnitude less accurate. 



12.1 from the IGM Lyman-a forest absorption 

The gas in the IGM is observed through H I Lya absorption in the spectra of 



all QSOs. Gunn & Peterson ||176|| discussed how redshift produces continuous 
absorption in the ultraviolet spectra of QSOs. Density fluctuations in the 
IGM trun this continuous absorption into the Lya forest absorption lines. 
The IGM fills the volume of space, and at redshifts z > 1 ||1 77fl it contains 
most of the baryons. 

The baryon density is estimated from the total amount of H I absorption, 
correcting for density fluctuations which change the ionization. The gas is 
photoionized, recombination times are faster in the denser gas, and hence 
this gas shows more H I absorption per unit gas. Using the observed ionizing 
radiation from QSOs, we have a lower limit on the ionizing flux, and hence a 
lower limit on the ionization of the gas. If the gas is more ionized than this, 
then we have underestimated the baryon density in the IGM. 

Three different groups obtained similar results ||178| , [|179|| , [ [18011 : flb> 



0.035/if 2 . This seems to be a secure lower limit, but not if the IGM is 
less ionized than assumed, because there is more neutral gas in high density 
regions, and these were missing from simulations which lack resolution. 

We do not have similar measurements at lower redshifts, because the 
space based data are not yet good enough, and the universe has expanded 
sufficiently that simulations are either too small in volume or lack resolution. 
Cen & Ostriker \ YF1 have shown that by today, structure formation may 



have heated most local baryons to temperatures of 10 5 — 10 7 K, which are 
extremely hard to detect [ 177| , [181 . 
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12.2 Clusters of Galaxies 



Clusters of galaxies provide an estimate of the baryon density because most 
of the gas which they contain is hot and hence visible. The baryons in gas 
were heated up to 8 keV through fast collisions as the clusters assembled. The 
mass of gas in a cluster can be estimated from the observed X-ray emission, 
or from the scattering of CMB photons in the Sunyaev-Zel'dovich (SZ) effect. 
Other baryons in stars, stellar remnants and cool gas contribute about 6% 
to the total baryon mass. 

The cosmological baryon density is obtained from the ratio of the baryonic 
mass to the total gravitating mass ||182| . Numerical simulations show that 



the value of this ratio in the clusters will be similar to the cosmological 
average, because the clusters are so large and massive, but slightly smaller, 
because shock heating makes baryons more extended than dark matter ||183| 



P-84JI . The total mass of a cluster, M t , can be estimated from the velocity 
dispersion of the galaxies, from the X-ray emission, or from the weak lensing 
of background galaxies. We then use Qb/ ^m— Mb/M t . The baryon fraction 
in clusters in the last factor is about 0.10/ifo 1 (SZ effect: | |185| |), or 0.05 — 
0.13/i TO 3/2 (X-ray: JTg§), or 0.1 lh^' 2 (X-ray: (TBJ, pg| ). Using Vt fi 



0.3 ±0.2 from a variety of methods ||189|| , we get fib— 0.03, with factor of two 



errors. These ilb estimates are lower limits, since there might be additional 
unobserved baryons. 

12.3 Local Dark Baryonic Matter 

The baryon density estimated in the Lya forest at z ~ 3 and in local clusters 
of galaxies are both similar to the that from SBBN using low D/H. This 
implies that there is little dark baryonic matter in the universe [ |190|| . This 



result seems conceptually secure, since there is little opportunity to remove 
baryons from the IGM at z < 3 or to hide them in dense objects without 



making stars which we would see [191], and the clusters are believed to be 



representative of the contents of the universe as a whole today. However, 
the numerical estimates involved are not yet accurate enough to rule out a 
significant density (e.g. 0.5 Qb) of baryonic MACHOS. 
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12.4 Simulations of the formation of Galaxies 



Ostriker (private communication) notes that the f2& can be constrained to a 
factor of two of that derived from SBBN using low D/H by the requirement 
that these baryons make galaxies. Semi-analytic models can also address 
the distribution of baryons in temperature and the total required to make 
observed structures (Frenk & Baugh, personal communication). 



12.5 CMB 

The baryon density can be obtained from the amplitude of the fluctuations 
on the sky of the temperature of the CMB. The baryons in the IGM at 
z ~ 1300 scattered the CMB photons. The amplitude of the fluctuations 
is a measure of fifji 2 , and other parameters. Published data favor large f^, 
with large errors, however dramatic improvements are imminent, and future 



constraints may approach or exceed the accuracy of from SBBN ||192 

EH- 



13 The Achievements of BBN 

Standard Big Bang Nucleosynthesis (SBBN) is a major success because the 
theory is well understood, close connections have developed between theory 
and observation, and observations are becoming more reliable. 

The early attempts to include physics in the mathematical model of the 
expanding universe lead to an understanding of the creation of the elements 
and the development of standard big bang theory, including the predictions 
of the CMB. 

The general success of SBBN is based on the robustness of the theory, 
and the resulting predictions of the abundances of the light nuclei. The 
abundances of 4 He, 7 Li and D can be explained with a single value for the 
free parameter rj, and the implied Qb agrees with other estimates. 

This agreement is used to limit physics beyond that in SBBN, including 
alternative theories of gravity, inhomogeneous baryon density, extra particles 
which were relativistic during BBN, and decays of particles after BBN. After 
decades owere f detailed study, no compelling major departures from SBBN 
have been found, and few departures are allowed. 
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Using SBBN predictions and measured abundances, we obtain the best 
estimates for the cosmological parameters r\ and Qb- 

The abundances of D, 4 He and 7 Li have all been measured in gas where 
there has been little stellar processing. In all three cases, the observed abun- 
dance are near to the primordial value remaining after SBBN. The D/H 
measured toward QSOs has the advantage of simplicity: D is not made after 
BBN, there are no known ways to destroy D in the QSO absorbers, and D/H 
can be extracted directly from the ultraviolet spectra, without corrections. 
There are now three cases of low D/H which seem secure. There remains 
the possibility that D/H is high in other absorbers seen towards other QSOs, 
but such high D must be very rare because no secure cases have been found, 
yet they should be an order of magnitude easier to find than the examples 
which show low D. 

We use low D/H as the best estimator of i] and the baryon density. SBBN 
then gives predictions of the abundance of the other light nuclei. These pre- 
dictions suggest that Y p is high, as suggested by Izotov, Thuan and collabo- 
rators. Low D also implies that 7 Li has been depleted by about a factor of 
two in the halo stars on the Spite plateau, which is more than some expect. 

The high f2b from SBBN plus low D/H is enough to account for about 
l/8th of the gravitating matter. Hence the remaining dark matter is not 
baryonic, a result which was established decades ago using SBBN and D/H 
in the ISM. 

The near coincidence in the mass densities of baryons and non-baryonic 
dark matter is perhaps explained if the dark matter is a supersymmetric 
neutralino |[L94|| . 

At redshifts z ~ 3 the baryons are present and observed in IGM with 
an abundance similar to f^. Hence there was no dark, or missing baryonic 
matter at that time. Today the same is true in clusters of galaxies. Outside 
clusters the baryons are mostly unseen, and they may be hard to observe if 
they have been heated to 10 5 - 10 7 K by structure formation. 

The number of free parameters in BBN has been decreasing over the years: 
Fermi & Terkovich gave nuclear reaction rates, the half-life of the neutron 
was measured, and then the number of families of neutrinos was measured. In 
standard BBN we are now left with one parameter, the baryon density, which 
is today measured with D/H using SBBN. When, in the next few years, this 
parameter is also measured, SBBN will have no free parameters. When free 
parameters can be adjusted to obtain consistency with the data, it is hard 
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to tell if a hypothesis is correct. The agreement between SBBN theory and 
measurement has grown stronger over the decades, as more parameters were 
constrained by independent measurements, and abundance measurements 
improved. This is the most convincing evidence that BBN happened and 
has been understood. 
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Figure 2: Abundances expected for the light nuclei 4 He, D, 3 He and 7 Li 
(top to bottom) calculated in standard BBN. New estimates of the nuclear 
cross-section errors from Buries et al. (1999a) and Nollet & Buries (1999) 
were used to estimate the 95% confidence intervals which are shown by the 
vertical widths of the abundance predictions. The horizontal scale, ry, is the 
one free parameter in the calculations. It is expressed in units of the baryon 
density or critical density for a Hubble constant of 65 kms _1 Mpc~ 1 . The 
95% confidence intervals for data, shown by the rectangles, are from Izotov 
and Thuan 1998a ( 4 He); Buries & Tytler 1998a (D); Gloeckler & Geiss 1996 
( 3 He); Bonifacio and Molaro 1997 ( 7 Li extended upwards by a factor of two 
to allow for possible depletion). 
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Figure 3: Optical spectrum of quasar 1937-1009, which shows the best example of 
primordial D/H. The top spectrum, from the Kast spectrograph on the 3-m telescope at 
Lick observatory, is of low spectral resolution, and high signal to noise. The continuum 
emission, from the accretion disk surrounding the black hole at the center of the quasar, 
is at about 6 flux units. The emission lines showing more flux (near 4950, 5820, 5940, 
6230, 6700 & 7420 A) arise in gas near the quasar. The absorptoin lines, showing less 
flux, nearly all arise in gas which is well separated from, and unrelated to the quasar. The 
numerous absorption lines at 4200 - 5800 A are H I Lya from the gas in the intergalactic 
medium. This region of the spectrun is called the Lya forest. This gas fills the volume 
of the intergalactic medium, and the absorption lines arise from small, factor of a few, 
fluctuations in the density of the gas on scales of a few hundred kpc. The Lya lines were 
all created by absorption of photons with wavelengths of 1216A. They appear at a range 
of observed wavelengths because they have different redshifts. Hence Lya absorption at 
5800A is near the QSO, while that at 5000 A is nearer to us. The abrupt drop in flux 
at 4180 A is caused by H I Lyman continuum absorption in the absorber at z = 3.572. 
Photons now at < 4180 A had more than 13.6 eV when they passed though the absorber, 
and they ionized its H I. The 1% residual flux in this Lyman continuum region has been 
measured in spectra of higher signal to noise (Buries & Tytler 1997) and gives the H I 
column density, expressed as H I atoms per cm -2 through the absorbing gas. The lower 
plot shows a portion of a spectrum with much higher resolution taken with the HIRES 
spectrograph on the Keck-1 telescope. We mark the Lya absorption lines of H I and D 
from the same gas. The column density ofQJ is measured from this spectrum. Dividing 
these two column densities we find D/H = 3.3 ± 0.3 x 10~ 5 (95% confidence), which is 
believed to be the primoridal value, and using SBBN predictions, this gives the most 
accurate measurements of rj and Qi,. 



